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Abstract. It has recently been established that imposing the condition of discrete 
holomorphicity on a lattice parafermionic observable leads to the critical Boltzmann 
weights in a number of lattice models. Remarkably, the solutions of these linear 
equations also solve the Yang-Baxter equations. We extend this analysis for the 
Zjy model by explicitly considering the condition of discrete holomorphicity on two 
and three adjacent rhombi. For two rhombi this leads to a quadratic equation in 
the Boltzmann weights and for three rhombi a cubic equation. The two-rhombus 
equation implies the inversion relations. The star-triangle relation follows from the 
three-rhombus equation. We also show that these weights are self-dual as a consequence 
of discrete holomorphicity. 



1. Introduction 

In some remarkable developments, a surprising connection has been uncovered between 
the notions of discrete holomorphicity and Yang-Baxter integrability [TJ [2j |3j HJ IS]- 
On the one hand, a parafermionic observable on the lattice is discretely holomorphic, 
i.e., the observable satisfies a version of the discrete Cauchy-Riemann equations. This 
requirement leads to a set of linear equations which can be solved to yield the Boltzmann 
weights of the model at criticality. The surprise is that, on the other hand, these are 
the same Boltzmann weights that follow by solving the star-triangle or Yang-Baxter 
equations. This connection has been observed by Rajabpour and Cardy [1] for the Zn 
model and by Ikhlef and Cardy [2] for the Potts model, the dilute 0(n) loop model and 
the dilute C% loop model. This list also includes the Ashkin- Teller model (HIS]- 

As yet there is no satisfactory explanation for this unexpected connection between 
discrete holomorphicity and Yang-Baxter integrability. Our aim in this paper is to 
explicitly establish that the star-triangle equation, and thus Yang-Baxter integrability, 
is a consequence of discrete holomorphicity. We will do this in the context of general 
rhombi on a Baxter lattice and the Z^ model. Our approach is thus algebraic. 
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The utility and power of the Yang-Baxter equation is well known. The importance 
and full power of discrete holomorphicity is becoming increasingly recognised. Setting 
up a discretely holomorphic observable has been a very useful step for the passage 
from discrete lattice models to the continuous functions of conformal field theory [6j[7]. 
Among other examples, discrete holomorphicity is a key ingredient in Duminil-Copin 
and Smirnov's \8\ remarkable and long sought after rigorous proof of the connective 
constant for self-avoiding walks on the honeycomb lattice. The exact value had been 
obtained 30 years ago by Nienhuis [9] in the n = limit of the 0(n) loop model on 
the honeycomb lattice. Nienhuis used various mappings between the vertex weights of 
different lattice models. The key point is that for n — all contributions from closed 
loops vanish, leaving only self-avoiding walks. The crucial first step in the arguments 
used by Duminil-Copin and Smirnov is the introduction of a parafermionic observable on 
the lattice which is discretely holomorphic. The second step involved the use of counting 
arguments in a domain of the hexagonal lattice. This step built on previous rigorous 
mathematical results obtained for decomposition of excursions across finite domains. A 
more formal proof has been given by Klazar [10] . The fact that such exact and rigorous 
results exist is no doubt related to the underlying integrability of the 0{n) loop model 
on the honeycomb lattice 

The arguments used by Duminil-Copin and Smirnov [8] were extended to the 0(n) 
loop model on the honeycomb lattice with a boundary by Beaton, Bousquet-Melou, 
de Gier, Duminil-Copin and Guttmann [12J to give a mathematically rigorous proof 
for the critical surface adsorption temperature of self-avoiding walks on the honeycomb 
lattice. This polymer adsorption transition corresponds to the n = limit of the 0(n) 
model at the special surface transition. The exact value for the adsorption transition 
had been obtained earlier by Batchelor and Yung [13] from the boundary Boltzmann 
weights obtained by solving the boundary version of the Yang-Baxter equation. Discrete 
holomorphicity at a boundary has been considered by Ikhlef [H] who was able to recover 
the boundary weights of the 0(n) loop model on the square lattice and to obtain a new 
set of boundary weights for the Zjv model. 

2. From discrete holomorphicity to the star-triangle equation 

For a Zjy symmetric iV-state spin model on a graph, the weight w(s a , Sb) of the 
interaction on an edge (a, b) is unchanged by the global transformations s r h-> ujs t 
and s r i — y s* of the spins. Here the spins take values from the iVth roots of unity, 
s r = u Qr , with q r G {0, 1, . . . , N— 1} and u = exp (2iri/N). For a review of these models 
and other important models contained therein, such as the Ising and Potts models, see, 
e.g., PS]. 

Since the weight depends only on the difference of the q r variables, we write it as 
W(q a — qb), or, by a slight abuse of notation, simply W(a — b). The Kramers- Wannier 
duality of the model [16] lets one define disorder variables fi? dual to the spins on the 
dual graph. The effect of inserting a disorder variable at the dual site f is to introduce 
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a string connecting the site to a site at infinity or at the boundary. Whenever the string 
intersects an edge of the original lattice, it modifies the weight on that edge by, for 
definiteness, lowering the spin q r to its right by one. 



p. .O* ■ ■ ■ ■ ■•' 
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Figure 1. A rhombic embedding of the covering lattice of a heterogeneous model 
showing a string attached to a disorder variable. The empty circles denote the spin 
variables, while the full circles denote the disorder variables. The location of the 
parafermion is marked with an x . 

We consider a rhombic embedding of the covering lattice, the union of the original 
lattice and its dual, onto the complex plane (figured]). Each elementary rhombus of this 
lattice contains an edge from the original lattice and an edge from the dual lattice. The 
parafermion introduced by Rajabpour and Cardy pQ 

1p( r> r) = exp ( — icr#( rj r)) • S r ■ [if 

lives on the midpoint of the edges. Here, the angle 9^) is the angle between the edge 
traversed from the spin to the disorder variable and an arbitrary but fixed axis. The 
real parameter a can be identified as the conformal spin of the observable. When the 
string enters a particular rhombus of opening angle a through the disorder variable Hx 
(figure [2]), imposing the condition of discrete holomorphicity 

^2 ^<r,f) &Z(r,r) = (1) 



on the contour sum of the parafermion around the rhombus in the counterclockwise 
direction we get, 

e Sg /iy e ^ ^ s a [iy ~\~ s a fj,<£ -\- e' e ' ^ Sg [i-^ = 

where the variables s a , [i X i s g-> are a ^ the corners of the contour, in that order, and 
Az is the difference between the complex coordinates of the two endpoints of the side 
of the rhombus on which the ip( r ,r) variable lives, traversed in the direction of traversal 
of the contour. 

As was explained in [TJ [14], we have taken the special case m = 1 of the more 
general definition of the parafermion 

Yfrfl = exp (-ia m 6( rt ?)) ■ S™ ■ ^ 

where the string attached to the disorder variable modifies the weight of an edge it 
crosses by lowering the spin to its right by m, and m £ {1, . . . , A^— 1}. These observables, 
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which are also discretely holomorphic at integrable critical points, are related to the one 
we consider by the global transformation s r i— > s™, and thus, it is enough to consider 
the m = 1 case only without loss of generality. 




, O ©-'. 



Figure 2. An elementary rhombus with opening angle a. The string enters the 
rhombus via the disorder variable The arbitrary axis is in the direction of the 
edge (a,X). 

With the introduction of 0(a) = e 1 ( 1 ~ a ) a ^ the condition takes the form 

0(-7r) s g fiy - 0(a) s a n Y + s a fix - 0(a - it) s g fi x = (2) 

Since, by our characterization of the disorder variables, 

fiy_ = W a {a-{g-l)) = W a {n a + 1) 
H% W a {a - g) W a (n a ) 

we have 

(0(-tt) - 0(a) uA) W a (n a + 1) + (u n " - 0(a - vr)) W a (n a ) = (3) 

which is a linear recurrence relation in the weights in terms of the bond varible n a . 

This relation specifies the weights up to a physically irrelevant constant. The 
symmetries, however, put constraints on the value of the unknown a (see, e.g., [T4"]). 
Since the weights are symmetric and periodic, we must have W a {n a ) = W a (N — n a ) for 
every n a , from which the recurrence relation gives 

(1 - 0(2a)) (uT 1 - 0(-2tt)) = 

Taking 

0(2tt) = u (4) 
makes the parafermion holomorphic for every opening angle a. Therefore, 

a = l-L+t (5) 



for any integer I. The case I = gives the Fateev-Zamolodchikov solutions jT7j. Since 
the subsequent calculations make use of relation (j3J) and not the actual value of a, the 
proof presented here is more general than merely proving that the Fateev-Zamolodchikov 
weights give us an integrable model. The angle a is proportional to (in fact, in [17], 
equal to) the spectral parameter and plays an analogous role. 

It should be noted that the weights are real because the ratio of the weights given 
by the recurrence relation ([3]) is real. 
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Figure 3. Rhombic embedding of an anisotropic square lattice. The edge with an 
opening angle a is horizontal, while that with an opening angle n — a is vertical in the 
original lattice. 



2.1. Self- duality 

We now show that the weights are self-dual. For this, we temporarily consider the 
embedding of an anisotropic square lattice (figure [3]). On the one hand, if the horizontal 
interaction edges are mapped onto rhombi of opening angle a, the vertical edges are then 
mapped onto rhombi with opening angle n — a. We thus have W a {n a ) = W n - a (n a ). On 
the other hand, taking the discrete Fourier transform of the recurrence relation ([3]) we 
have 

(0(-tt) - 0(-tt - a) ou k " +1 ) W a (k a + 1) + {u h « - cf>(-a)) W a (k a ) = 

Using (j3J) to replace the extra u in the first term, we see that the dual weights satisfy the 
same equation with a replaced by n — a. Since our unitary Fourier transform preserves 
the norm 

E t^K)] 2 = 1 (6) 

ria&N 

this proves 

W a (n a ) = W n . a (n a ) = W a (n a ) (7) 
and captures the crossing symmetry of the problem. 

2.2. Inversion relations 

Now we consider another rhombus with angle which sits on top of the rhombus 
considered so far (figure H]) . The new rhombus comes with spin Sb and disorder variable 
fig and shares the edge (g,Y) with the previous one. We multiply its holomorphicity 
equation analogous to (T5]) by </>(—/?) to orient it correctly and then add the two equations 
to get 

- <f)(a) S a fly + SaUx - <K a - S g fix 

+ 0(-vr - 0) SgHz - s b fiz + 0(-/3) S b [ly = (8) 
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Figure 4. Two adjacent rhombi. 

The contributions from the common edge cancel as it is traversed in opposite directions 
for the two sums. Using 

fj.2 _ W p (b-(g-l)) _ W p {n h +l) 
Hy Wp(b-g) W p (n h ) 

the equation (jSJ) gives a quadratic relation 

( w «-0(a-7r) u 9 )W a (n a )Wp(n b ) 

+ co b - 0(a) u a ) W a (n a + 1) Wp(n b ) 

+ (0(-tt - 0) u 9 - oo b ) W a (n a + 1) W p (n b + 1) = (9) 

in the weights. 

This relation contains the inversion relations (see, e.g., [18]) for the model. Put 
a = b, so that n a = n b = n, and (3 = —a, then the equation reduces to 

W a {n) W^ a {n) = W a (n + 1) W. a {n + 1) 
thus, the product is independent of n, 

W a (n)W. a (n)=g(a)g(-a) (10) 
The other relevant inversion relation 

w *+o( a - 9) W„-e(b -g) = p{6) 5 a , b (11) 

follows from this relation and self-duality by expanding the left hand side in Fourier 
coefficients and gives 

p(6) = N g(6) g(-6) 

for our normalization. 

2.3. Star-triangle relation 

We proceed to add the third rhombus with angle 7 = 2-k — (a + 0) to make the star 
(figure [5]). Two of the edges, (g, Z) and (g,X), are shared and a new spin variable s c is 
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Figure 5. Three adjacent rhombi making the star. 

introduced. Multiplying the holomorphicity equation analogous to for this rhombus 
by 0(27r — a) to orient it correctly, we add it to the contour sum OH]) to yield 

(p(a — 3n) Sgfix — (f>(a — n) s g \x x 

+ S a Hx ~ 0(a) S a 

+ 0(-/3) S b yUy - Sfc/i^ 

+ 0(a-27r)s c ^-0(-/9)a c ^ = O (12) 

Some features of this equation warrant attention. The disorder variable \x x is not the 
same as p, x , as might be naively expected. A possible interpretation of this is that the 
string has gone around the spin s g once, and thus the parafermion has acquired a phase 
factor 0(27r), corresponding to the lowering of the value of g, that is, q g , by one in the 
configuration. Hence, the contributions from the edge (g, X) do not cancel. 
Making use of 

U x _ W 7 ( c -(3-l)) _ W 7 (n c + 1) 

fly W 7 (c-g) W 7 (n c ) 

we have a cubic equation in weights which we omit for clarity. However, summing over 
g, we arrive at 

{u a -<P(-f3)uj c ) £ a ,p(a,b,c) 

+ - 0(a) a/ 1 ) £ a)j9 (a + l,&,c) 

+ (0(a - 2tt) cu c - u b ) £ a! p(a + 1, b + 1, c) = (13) 

Here 

£« )j9 (a, b, c) = W ^ - 9) Wpib - g) W 7 (c - a) (14) 

is the left hand side of the star-triangle relation. The terms involving s g cancel because 
of the constraint @ on a, and therefore, the expectation value of the parafermion is 
still independent of the path, as is that for the disorder variable fx x , and consequently 
the contour sum around the star still vanishes. 

Remarkably, considering the contour sum around the triangle (figure , we get 

(w'-^w 6 ) 9W(a,6,c) 
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+ (<j>(-P) co b - 0(a) tu a ) JH«^(o + 1, b, c) 

+ (0(a - 2tt) w c - c/) *RaA a + 1, & + 1, c) = (15) 

where 

^(o, b, c) = W^(a - b) W„- a (b - c) W„- P {c - a) (16) 

is proportional to the right hand side of the star-triangle relation. The two sides of the 
equation therefore obey the same recurrence relation. 

To show that the star-triangle relation follows from these relations, we have to 
prove that the two solutions are linearly dependent. To this end, we consider a general 
solution of the recurrence relation f a; p(a,b,c). Taking the complex conjugate of the 
relation and noting both f a ,p(a,, b, c) and a are real, we multiply the resulting equation 
by (j){a - (3) u a+b . The result, 

(<f>(a -(3)u b - <P(a) u a+b - c ) f a> p(a, b, c) 
+ {<P{a)uj a -<P{-f3)uj b ) f a ,p(a + l,b,c) 

+ (0(2tt - (3) u a+b ~ c - 4>(a - 13) u a ) f a ,p(a + 1, b + 1, c) = (17) 

when added to the recurrence relation, cancels the middle term and gives the ratio 

fa,p(a + l,b+ l,c) = f a ,p(a,b,c- 1) ^ ig ^ 
f a! p{a,b,c) f ai p(a,b,c) 
independent of the function f a ^. The shift in the arguments is justified by the Zn 
symmetries. Since this ratio is the same for both Dl a ,p, and £ a ,p, 

£ ai p(a,b, c) _ £ ai p(a,b,c- 1) 
fRa^(a,b,c) di a> p(a,b,c-l) 

The ratio of the two sides is thus independent of the spin c. Similar elimination of the 
other two terms in the recurrence relation shows that the ratio is independent of a and 
b as well, as is demanded by symmetry^ We therefore have the star-triangle relation 

& at fi (a, b, c) = R a; p JR^ (a, b, c) (20) 

| Jacques Perk informed us that this style of argument has been used to prove that the Boltzmann 
weights of the more general chiral Potts model satisfy the star-triangle relation [19] . 
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where R a ^ depends on the spectral variables only. 

As is well-known, the inversion relations f lT0|) - f|TT|) and the star-triangle relation 
( 1201 together are sufficient to establish the existence of commuting transfer matrices 
parametrized by, in our case, the angle of the embedded rhombi, or in other words, 
Yang-Baxter integrability. By the embedding onto the complex plane, the schematic 
diagram of the Yang-Baxter equations acquires a geometric meaning of rearrangements 
of the elementary rhombi. Our analysis shows that the contour sum around a domain 
picks up only factors independent of the configuration under such rearrangements. 



2.4- Homogeneous lattices 
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Figure 7. Rhombic embedding of the three homogeneous lattices. 



It is interesting to note that the construction presented above is local in that it does 
not depend on many of the properties, e.g., the coordination number, of the underlying 
graph. To see the consequences of this lattice independence, we consider the three 
archetypical homogeneous lattices, the square, the triangular and the honeycomb (figure 
[7]). The embedding of these lattices onto the complex plane tiles the plane with a being 
7r/2, 7r/3 and 27r/3 respectively. Evaluating the recurrence relation (|3]) at these points 
gives an instant derivation of the known critical weights. For example, for the Ising 
model, we have for W(1)/W(0), the known values \/2 — 1, and 2 — \f?> for the 

three lattices respectively [20J. 



3. Concluding remarks 

In this paper, we considered the implication of the condition of discrete holomorphicity 
on two and three adjacent rhombi in the context of the lattice Z^r model. For two 
rhombi this led to the quadratic equation (jSJ) in the Boltzmann weights. This equation 
was shown to imply the known inversion relations ffTUj) and (fTTj) for this model. Note 
that Cardy [21] has effectively established the inversion relations in the u — v — > 
limit of the Yang-Baxter equation. For three rhombi we obtained the cubic equation 
(JT2l) in the Boltzmann weights. In establishing this equation the lattice parafermion 
picks up a crucial phase factor, with the expectation value of the parafermion still 
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independent of the path. The importance of such a phase factor has been highlighted 
in the topological context of the loop models by Fendley [22J. Here we have shown 
that the star-triangle relation ( 1201) follows from the three- rhombus equation ( Tl2|) . In the 
discrete holomorphic approach the two-rhombus equation (jSj) and the three-rhombus 
equation ( TT2"]) can thus be considered as analogues of the two- and three-body conditions 
for integrability. However, the simplicity of the discrete holomorphic approach is that 
ultimately these conditions are equivalent to the one- rhombus equation (j2]). Indeed, 
one can push this argument further by building up a transfer matrix from the rhombi 
and showing that commuting transfer matrices can be established as a consequence of 
discrete holomorphicity, bypassing the use of the Yang-Baxter equation. Our results 
lend further impetus for using discrete holomorphicity as a tool for investigating lattice 
models at critical points. 
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